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ABSTRACT 



Aims. We create a model for recovering the intrinsic, absorption-corrected surface brightness distribution of a galaxy and apply the 
model to the nearby galaxy M 3 1 . 

Methods. We constructed a galactic model as a superposition of axially symmetric stellar components and a dust disc to analyse the 
intrinsic absorption effects. Dust column density is assumed to be proportional to the far-infrared flux of the galaxy. Along each line 
of sight, the observed far-infrared spectral energy distribution was approximated with modified black body functions corresponding 
to dust components with different temperatures, thereby allowing us to determine the temperatures and relative column densities of 
the dust components. We applied the model to the nearby galaxy M 31 using the Spitzer Space Telescope far-infrared observations for 
mapping dust distribution and temperature. A warm and a cold dust component were distinguished. 

Results. The temperature of the warm dust in M 31 varies between 56 and 60 K and is highest in the spiral arms, while the temper- 
ature of the cold component is mostly 15-19 K and rises up to about 25 K at the centre of the galaxy. The intensity- weighted mean 
temperature of the dust decreases from T ~ 32 K in the centre to T ~ 20 K at ~ 7 kpc and outwards. The scalelength of the dust disc 
is (ao)dust ~ 1.8 (ao)stais- We also calculated the intrinsic U, B, V, R, I, and L surface brightness distributions and the spatial luminosity 
distribution. The intrinsic dust extinction in the V-colour rises from 0.25"' at the centre to 0.4"' - 0.5"' ^tR^6 - 13 kpc and decreases 
smoothly thereafter. The calculated total extinction-corrected luminosity of M31 is = (3.64 ± 0.15) • IO'^Lq, corresponding to 
an absolute luminosity Mb = -20.89 ± 0.04 mag. Of the total 6-luminosity, 20% (0.24 mag) is obscured from us by the dust inside 
M 3 1 . The intrinsic shape of the bulge is slightly prolate in our best-fit model. 
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The importance of correcting astronomical observations for the 
impact of intervening dust has been known for a long time. 
Several direct and indirect methods have been developed for es- 
timating light absorption and reddening by interstellar matter. 
However, because of the complexity of this process, the influ- 
ence of dust is often ignored in studies of individual galaxies, 
^ as well as extensive statistical samples. Unfortunately, ignor- 
ing dust effects seriously undermines the reliability of the re- 
sults, since it may lead to substantially wrong estimates of the 
shapes, luminosities , colours, and masses of the stellar com- 
ponents of galaxies dPriver et alj|2007l: iPadilla & StraussllIOOSt 
iGraham & Worlevll2008h . 

Two significantly different approaches can be used to es- 
timate intrinsic dust effects in external galaxies. Firstly, light 
extinction and reddening can be traced directly by comparing 
the properties of affected and unaffected fields within or be- 
hind a given galaxy. For example, in the presence of a partly 
occulted background galaxy, the optical thickness and extinc- 
tion cur ve of dust i nside the foreground galaxy can be measured 
directly dHolwerda e t al. 2007). Unfortunately, suitable pairs of 
galaxies are very rare. Alternatively, extinction can be estimated 
by counting background galaxies seen through the foreground 
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galaxy and by comparing their abundance to unobscured fields 
(Gonzalez et al. 1998, 2003; Holwerda et al. 2005). The statisti- 
cal nature and difficulties in detecting galaxies behind brighter 
regions (e.g. the bulge) of the foreground galaxy tend to reduce 
the applicability of this method. 

The second approach exploits the law of conservation of 
energy - optical and ultra-violet radiation absorbed by dust 
has to be reradiated at the infrared (IR) and longer wave- 
lengths. Data collected by the Infrared Astronomical Satellite 
(IRAS), the Infrared Space Observatory (ISO), the Spitzer Space 
Telescope, and the Submillimetre Common-User Bolometer 
Array (SCUBA) have significantly improved our understanding 
of the spatial distribution, temperature, and physical properties 
of interstellar dust grains, enabling the construction of radiation 
transfer models even down to the level of individual dust grains. 
In addition to revealing new facts about cosmic dust itself, these 
models enable restoration of the intrinsic photometric properties 
of galaxies for which direct methods of estimating light extinc- 
tion are not applicable. 

Since the thermal radiation from dust strongly depends on 
dust temperature, far-IR emission maps taken at several wave- 
lengths are required for calculating dust density and for estimat- 
ing the extinction of starlight. The first proof that dust tempera- 
ture may significantly vary within a galaxy was noted from the 
IRAS data by Helou (1986) and Lonsdale & Helou (1987). It 
was suggested that, in addition to the warmer dust associated 
with star-forming regions, a cooler and more diffuse dust com- 
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ponent exists. Higher angular resolution far-IR maps of several 
galaxies by the ISO satellite indicated that the spatial distribu- 
tions of t he cold dust and the warm dust are different (Haas e t al 
1998; T uffs & Gabrielll2003HEngelbracht et al.ll2004tlHinz et al 
2004) . It was also discovered that the cold dust disc is us ually 



more extended than the stellar disc (iPopescu & Tuffsll2003h 

By approximating the IRAS and SCUBA measurements 
of dust spectral energy distributions in 18 galaxies with two- 
component dust models, IVlahakis et al.l (|2005) found that cold 
dust has temperatures 17-24 K and warm dust 28-59 K. By us- 
ing the Spit zer, SCUBA and IRA S observations for a sample of 
10 galaxies IWillmer et al.l (|2009) derived the temperature range 
for the cold dust 1 8-24 K and for the warm dust 52-58 K. The ra- 
tio of the cold dust mass to the warm dust mass in these galaxies 
is typically 100-2000. These values were derived by averaging 
across each of the galaxies; dust temperature may be consider- 
ably different in specific regions. 

To determine the total mass of du st in g alaxies, Li & Drainp 
(1200 lb . IWeingartner & Draind (1200 Ih . and iDraine & Lil (l2007h 
proposed a model considering three kinds of dust grains and 
the radiation field, in which dust temperature variations from 
galaxy to galaxy are regulated by the effective intensity of the 
radiation fi eld. This model has been successfully applied to 17 
galaxi es by Draine et al. j 2007h and more recently to 57 g alax- 
ies by 'Munoz-Ma teos et all (l2009l) . iMontalto et"al] (120091) have 
used this model for studying M 3 1 . 

More sophisticated analyses, yielding dust temperature as 
a function of position in a galaxy and thus enabling a de- 
tailed comparison to far-IR maps are carried out within radia- 
tive t ransfer models of stellar radia t ion in a dusty environment 
(e.g. ISautv eFai] 119981; ISilva et all 119981; iBianchi et all 120001: 
iPopescu et al. 20001; Tuffs et al.1 2004 ) Spa tial resolution can 
reach 20 pc in these models "(iBianchil 120081) . However, this is 
still not enough to resolve individual dust clumps. The problem 
can be solved by separating the clumpy warm dust compo nent 
from the diffiise com ponent, as done by jSilva et al.l (1 19981) and 
IPopescu et all (|2000|) . 

In this paper, we create a three-dimensional galaxy model 
with axi-symmetric stellar populations. To estimate extinction 
and to restore the intrinsic luminosity and colour distributions 
of a galaxy, we add a dust disc to the model. We apply the 
model to the nearby spiral galaxy M31, using the Spitzer far- 
IR maps (supplemented by the IRAS observations) to determine 
the dust distribution. The distribution of the dust temperature is 
estimated by approximating the far-IR spectral energy distribu- 
tion with modified black body functions assuming the presence 
of a colder and a warmer dust component. 

The following general parameters of M31 have been ap- 
pUed in our calculations: th e inclination angle has been take n 
77.5° (IWalterbos & Kennicut t 1988; de Vaucouleurs et al.lll99ll); 
the major axis position angl e is 38.1°(lWalterbos & Kennicura 

the distance to M31 is 
corresponding to the scale 



119871; Ferguson et al] |2002|) and 



I2005h . 



785 kpc (Mc Connachie et al. 
1' =228pc. 



2. Construction of the model 

The physical properties of dust and the extinction of the stellar 
light of a galaxy are determined by an interplay between the spa- 
tial radiation field and the dust grains at each location within the 
galaxy. Calculations of the intrinsic extinction of a galaxy should 
thus be based on the spatial luminosity distribution of the galaxy 
and the spatial density distribution of the dust grains. 



2. 1 . Density distribution model for the stellar components 

In this subsection, we develop a sufficiently flexible model for 
describing the spatial distribution of the luminosity of a galaxy. 
A two-dimensional projection of the model can be compared to 
the observed surface brightness distribution and model parame- 
ters can be adjusted. 

The model galaxy is given as a superposition of its individual 
stellar components. Each visible component is approximated by 
an ellipsoid of rotational symmetry with a constant axial ratio q; 
its spatial density distribution follows the law 



1(a) = 1(0) exp 



a 



I IN 



(1) 



where 1(0) - hLKAnqa^^ is the central density and L is the com- 
ponent luminosity; a - -\jr^ + z^/q^, where r and z are two cylin- 
drical coordinates; qq is the harmonic mean radius that charac- 
terises the real extent of the component, independently of the 
structure parameter A^. The coefficients h and k are normalis- 
ing parameters, dependent on A^. The definition of the normal- 
ising paramet ers h and k and the ir calculation is described in 
appendix B of iTenies et aL ( 1994i). The luminosity densi ty dis- 
tribution ^ proposed by lEinastoTd 19651) is similar to the lSersij 
( 1196 8) law for surface densities. Differences between the Sersic 
law and the two - dimen sional projection of Eq. [T]can be seen in 
iTamm & Teiiiesl (|2005|) . 

Observations have demonstrated that the stellar disc can have 
a toroidal form in some galaxies, i.e. it does not continue to 
the centre. As suggested by Einasto (' Einastoll969l ; lEinasto et al] 
il98Q), the spatial density of a disc with a central hole can be 
expressed as a sum of two spheroidal mass distributions 



'di.sc(a) = l+(a) + l-(a). 



(2) 



both of which can be approximated with the exponential law 
([T]i. Adopting a model disc with zero density at r = and non- 
negative density at Zdisc(fl) > gives the following relations be- 
tween the parameters of the components 1+ and L: aq- = Kao+, 
L - -K^L+, q^ - q+lK, where k < 1 is a parameter that de- 
termines the relative size of the hole in the centre of the disc. 
Structural parameters and should be identical. 

Density distributions of all visible components are projected 
along the line of sight and their sum yields the surface brightness 
distribution of the model 



L(A) = 2 



oo 



lj(a)a da 
(a^-A^yi^' 



(3) 



where A is the major semi-axis of the equidensity ellipse of the 
projected light distribution and Qj are their apparent axial ratios 
2^ — cos^ i + q^ sin^ The inclination angle between the plane 
of the galaxy and the plane of the sky is denoted by /. The sum- 
mation index j designates each visible component. 

To approximate the luminosity distribution of a given galaxy, 
an arbitrary number of components can be used. Obviously, a 
larger number of components enables a better approximation of 
observations; however, to keep the model physically justified, 
additional information (e.g. observations of metallicity and kine- 
matics distribution) should be used for distinguishing different 
components and for determining their number and extent. For 
the majority of nearby regular galaxies, three to five stellar com- 
ponents with density distribution ([T]) should be enough for de- 
scribing the luminosity distribution. 
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Fig. 1. Geometry of the dust disc model. The line-of-sight optical 
depth for three positions A, B, and C is indicated. The figure is 
illustrative only. 



The model parameters of each component are determined by 
a least-squares approximation of the observed surface brightness 
distributions with the model profiles via Eq. [3] Uncertainties 
of the model can be estimated by using the partial sec- 
ond derivatives of the sum of the least-square differences 
jBevington & Robinsonll2003l) . 



2.2. The dust disc component 

Let us now implant a dust disc component into our galaxy model, 
allowing us to calculate absorption along each line of sight 
through the galaxy. 

Equation [3] will give the luminosity density along a line of 
sight, if no dust extinction is included. The intrinsic luminosity 
density for each component can be calculated from the equation 



2 [l(r,Zj)e-''-'''^^] 



L(X, Y) 



Zl,2 



•J sm; 



-rdr. 



tan I 



(4) 



(5) 



where /(r, z) denotes the spatial luminosity density ([T]l of the 
galaxy component, L(X, Y) is the corresponding surface bright- 
ness distribution, X and Y are coordinates in the plane of the sky, 
and T{r, z) is the optical depth along the line of sight (see Fig.[D- 

Optical depth T{r,z) is zero between the observer and the 
dust disc and T(r,z) - T^ax behind the dust disc. T^ax is the total 
optical depth along a given line of sight and thus a function of 
(X, Y). Inside the dust disc, T(r, z) is between and Tmax- 

The optical depth between the observer and each point inside 
the dust disc can be written as 



T(r,z)^T^AX,Y) 



J ndust(i) ds 

A 

"c ' 

/ «dust(«) ds 

A 



(6) 



where integration is done along the line of sight s and rtdusii-^) is 
the spatial number density of the dust disc. 

To gain a simple but nevertheless flexible model for the dust 
disc, we adopt the following form for ndust^ 



«dust('", z) = fdu^r) ■ n(a). 



(7) 



where n{a) is the exponential law ([TJ and fdusM describes devi- 
ations from the exponential law ([T]i. The spatial number density 
of dust decreases exponentially in the z-direction according to 



Eq. [T] Multiplication by fdustir) allows density behaviour to be 
more flexible along the r-direction. 

Assuming that extinction is proportional to the dust column 
density and the latter is proportional to the far-IR flux, the pa- 
rameters of n(fl), the shape of /du.st('") and subsequently the map 
of Tmax(-'^, Y) can be derived. The relation between reddening and 
the far-IR flux can be expressed as 



E(B - y) = pD, 



(8) 



where D is a measure of the dust column density and p is 
a calibration constant dSchlegel et alj [19981) . Expressing D as 
the 100//m flux of dust with a reference temperature 18.2 K, 
ISchlegel et aPd 19981) have estimated p = 0.0 16+0.004 by study- 
ing the colours of bright elliptical galaxies as seen through the 
Galactic dust. 

To make Eq. [8] applicable for other galaxies, we need to take 
in account the difference between the lOOyum flux of dust at the 
temperature 18.2 K and the actual flux of dust with a temperature 
r at a wavelength A. According to the modified Planck function 



D(X, Y) = FAX, Y) 



{lOOiumfBAA, T(X,Y)y 



(9) 



where F^X, Y) is the measured far-IR flux map at a wave- 
length A and T(X, Y) is the temperature map, is the black 
body function and /? is the dust emissivity index. According to 
Vlahakis et al. (2005) we take /3 = 2.0. 

Using Eq.|8]and the relation Ay = /? -£'(B-y) and converting 
the derived extinction into optical depth, we reach a map of the 
total optical depth t^^^^X, Y) 

In(lO) 

Tmax,i(^, Y) = -^RpD(X, Y). (10) 

Z.J Ay 

IXilouris etaf] (11999 ^ have found that in several galaxies, the ex- 
tinction law agrees well with that of the Milky Way. This is also 
the case for M31 as i ndicated by the colours of its globular 
clusters (iBarmbv et al.|[2000.) . We can thus use p - 0.016 and 
R - 3.1; the extinction ratio A^/Ay can be determined from the 
Galactic extinction law. 

Equation[TO]gives the total extinction map in the plane of the 
sky. To calculate the actual intrinsic extinction, we also need the 
spatial density distribution of the dust disc ndu.st('", z) for Eq. |6] 
We can compare the projection of the model space density dis- 
tribution (calculated according to Eq. [3]) to the observed dust 
column density distribution (i.e. the far-IR maps) to determine 
the function /dust('") and the parameters for n{d). The observed 
column density of the dust disc is thereby converted into its spa- 
tial number density. More precisely, we can only determine the 
shape of the number density distribution of the dust disc; its ab- 
solute calibration n(0) remains unknown. In the present case, the 
calibration constant is cancelled out in Eq.|6]and is not required 
for our model. 

In the final step, the derived optical depth map Tmax,^(-'^, Y) 
and the function ndust('", z) are inserted into Eq. |6] and the 
extinction-corrected surface brightness along each line of sight 
can be calculated from Eq.|4l 

The relation between Tmax,^(^, Y) and reddening in Eq.fTOlis 
calibrated according to the reddening of extragalactic sources by 
the dust of the Milky Way, thus including both extinction and 
scattering. On the other hand, the luminosities calculated from 
Eq.lH include dust attenuation, but ignore the scattering of pho- 
tons, leading to an overestimation of stellar emission in regions 
with high dust density. However, we consider the additional un- 
certainties caused by scattering to be lower than other potential 
error sources in r determination. 
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Table 1. Used photometrical data. 



Reference 


Filter Remarks 




system 



de Vaucouleurs (1958 
Hoessel & Melnick (1980 



Hodge & Kennicutt (1982) 
Kent (1987) 

Walterbos & Kennicutt (1988) 
Pritchet & van den Bergh (1994) 




Guhathakurta et al. (2005) 
Irwin et al. (2005) 

Barmbv et al. (2006) 



VI 



ellipses 
major and 
minor axis 
ellipses 
ellipses 
ellipses 
star counts 
(minor axis) 
star counts 
(minor axis) 
ellipses and 
star counts 
(minor axis) 
ellipses 



Note: 
files. 



'ellipses' denote elliptically averaged surface brightness pro- 



3. Applying the model to M31 

3.1. Optical data 

A large collection of surface brightness distribution measure- 
ments is available for modelling the s tellar luminosity distribu- 
tion in M 3 1 ; see iTenjes et~an d 19941) for references to earlier 
data. For the present work we have used the measurements of 
U, B, V, R, and / distributions excluding some of the oldest, 
mainly photometrical data; however, for a more complete cover- 
age of B surface brightness distribut ion, we have also inclu ded 
the photoelectrical measurements bv .de VaucouleursI (Il958l) . As 
an example of nearly absorption-free observations, we have 
added the 3.6ium me asurements by the Spitzer Space Telescope 
dBarmbv et al.l200^ to our dataset, which roughly correspond to 
the Johnson-Cousins L filter Since we did not intend to model 
the nuclear part of M31, data at r < 0.3 kpc were neglected. 
Table [1] presents references to the observations and the corre- 
sponding filter systems used in this paper Some of the data are 
one-dimensional measurements along either or both of the axes 
of the galaxy, others are elliptically averaged contours of the 
galaxy. Measurements along the axes are useful for an accurate 
determination of the three-dimensional shape of the stellar pop- 
ulations, while ellipses provide the average surface brightness 
at each distance from the centre. Of the minor axis data, only 
those from the side closer to us (the southeast direction) have 
been used and are presented on figures throughout this paper. 
Different R and / colour system profiles were transfe rred into the 
Johnso n-Cousins system, using the calibrations by Frei & GunnI 
(|1994|) . All measurements wer e corrected for extinction by dust 
of the Milky Way according to lSchlegel et all fl998 ). 

3.2. Infrared data 

Our analysis of the far-IR emission from M 3 1 is based on 
the highest spatial resolution mappings currently available, the 
Spitzer Space Telescope MI PS camera imaging at 24, 70, and 
160/zm dGordon et allboO^ . supplemented by the IRAS satel- 
lite detections at 100/zm. The only measurements available at 
longer wavelengths are those by the DIRBE instrument aboard 
the COBE satellite at up to 240 jum. Because of their very low 
angular resolution, the DIRBE observations were not suitable 
for this work. 
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Fig. 2. Examples of approximations of the observed Spitzer 
(filled circles) and IRAS (empty squares) data with the modified 
black body functions at three random locations: (a) in the bulge 
region; (b) in an interarm region; (c) in a spiral arm. In each 
panel, the black dashed lines represent the warm and the cold 
component and the red solid line is the sum of the two modified 
black body functions together with its uncertainty region. 



The Spitzer images were downloaded using the Leopard 
toolkit developed at the Spitzer Science Center. An inspection 
of the pipeline-created (post-BCD) mosaic images revealed un- 
even background levels and the presence of a number of in- 
strumental defects, making them unusable for our purposes. We 
started from the basic-calibrated data (BCD) and conducted the 
necessary corrections, stacking and mosaicing using the stan- 
dard procedures within the IRAF and MOPEX software. The 
first three frames of each 24 yccsx observing run, suffering from a 
shorter exposure time, were discarded. Depending on the fields, 
either a constant or a low order polynomial background was sub- 
tracted from the images. The IRAS imaging was retrieved from 
the NASA/IPAC Infrared Science Archive. The resolution of the 
IRAS data are only 90 arcsec/pixel; they were not used during 
the final approximation of the spectral energy distribution. The 
final images were sampled to 8 arcsec/pixel, which is the sam- 
pling of the 160 yum observations; this also became the spatial 
resolution of our dust model. 

Uncertainties were estimated considering three principal er- 
ror sources: the characteristics of the telescope detectors, the pre- 
cision of the absolute calibration and the uncertainty of the back- 
ground level. An estimate for the response uncertainty of each 
detector pixel is indicated by the corresponding standard devi- 
ation maps; these were stacked and resampled by the MOPEX 
software during the mosaicing process. The conservative esti- 
mates of the uncertainti es of the instrumental flux calibration of 
Spitzer were taken from Engelbracht et al.l ( l2007h . lOordon et aTl 
(2007), and Stansberry et al (2007 ): 4%, 7% and 12% for 24, 
70, and 160//m, respectively. Bei chman et alJdl988[) give 2 % as 
the maximal variation in the IRAS absolute calibration. Possible 



E. Tempel et al.: Dust-corrected photometry of M31 



5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 





I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



I I I I I I III I I I I I I I I I I I I I I I I I I I I 



I I I I I I 







g 


61 




4 


S3 




1 = 






^. 


57 




' -2 


55 




-4 


2 










4 






■0- £ 


1 




it) 






-2 






-4 







-€ 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I - I I I I I I I 



I I I I I I I I I III I I I I I I I I I I I I I I I I I I I I I 



- t 

I 



0.5 ^..^ 



0.1 



■15 -10 -5 5 10 15 
X (kpc) 



Fig. 3. Temperature and intensity maps of the warm dust com- 
ponent. In panels (a) and (b) the warm dust temperature and its 
Icr-errors are given. In panels (c) and (d) the intensity of the 
warm component at 100 yum and its Icr-errors are given. 
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Fig. 4. Temperature and intensity maps of the cold dust compo- 
nent. In panels (a) and (b) the cold dust temperature and its Icr- 
errors are given. In panels (c) and (d) the intensity of the cold 
component at 100 yum and its Icr-errors are given. 



background errors were estimated from deviations between the 
background flux density levels of neighbouring data stripes. 
Combining these three error sources, uncertainty maps were pro- 
duced for each image, to be later used for calculating possible 
errors in our extinction analysis. In addition, the uncertainty of 
disentangling the contributions of M 3 1 and the Milky Way from 
the far-IR maps has to be kept in mind. According to the IRAS 
maps, the sky bri ghtness at lOO/im ar ound M 3 1 varies between 
2 and 5MJysr ' jSchlegel et al.llT998i) . This also gives an idea 
about a possible error resulting from wrong addressing of dust 
in the direction of M 3 1 . We have not included this foreground 
uncertainty in our error estimates. 

To probe the repeatability of our image processing and pho- 
tometry, let us compare our estimate for the integral IR flux to 
some earlier measurements. The total sky-subtracted flux within 
the outermost ellipse fitted to the galaxy with the IRAF/STSDAS 
task ellipse is 119+7Jy, 1200+1 10 Jy, and 7800±1000Jy for 
the 24, 70, and 160 /zm images, respectively. These corre- 
spond rather wefl to the value s 107+lOJy, 9 40±188Jy, and 
7900±1580Jy measured by Gor don et al] (|2006) inside a rect- 
angular aperture (2.75°x 0.75°). The variations probably come 
from different methods used for background subtraction and 
from the aperture selection. 

3.3. Mapping dust temperature and optical deptin 

Following Eqs. |9] and (TO] we need a far-IR emission map and 
a map of the dust temperature. Since we are relying on the rela- 
tio n between reddening a nd the far-IR flux at 100 ;um determined 
bv lSchlegel etal.1 (1 19981) for the Milky Way, the far-IR map has 



to correspond to the lOOyum flux. Although the IRAS satellite 
has observed at this wavelength, we derived the necessary map 
by interpolating the Spitzer observations because of their higher 
spatial resolution. 

iGordon et al.l (l2006l) have shown that a two-component dust 
model gives a good fit to the far-IR spectral energy distribution 
of M31. The warmer (T ^ 59 K) dust component is about 10^ 
times less massive than the colder (T ^ 17K) component, but it 
might still cause significant light extinction in specific regions, 
since it is considerably more clumpy and more concentrated in 
the galactic centre and in the spiral arms. Therefore, we have fit- 
ted the far-IR emission within each line of sight with two mod- 
ified Planck functions, representing these two dominant temper- 
ature components. 

To determine the intensity and temperature of both of the 
dust components, the modified Planck functions were fitted to 
the 4-bin spectral energy distribution observations (24, 70, 100, 
and 160 ^m) in three stages for each line of sight (i.e. within each 
model pixel). The main contribution to the warm component 
comes from the star-forming regions. We assumed that the tem- 
perature of dust in these regions varies only slightl y and fixed the 
temperature of the warm component at T = 59 K (iGordon et al.l 
2006) during the first fitting stage; the other parameters were de- 
rived using the Spitzer data. During the second stage, the IRAS 
data were added and the temperature of the warm dust was ad- 
justed at lower spatial resolution. Finally, IRAS data were omit- 
ted again and the intensities of the dust components and the 
temperature of the cold dust were refined at higher spatial reso- 
lution. The uncertainties of the derived intensities and tempera- 
tures were estimated by using the Monte Carlo method, allowing 
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Fig. 5. Four-folded optical depth map of M 3 1 , derived from the 
cold dust component using Eq. \W\ Panel (a) - the observed op- 
tical depth; panel (b) - the modelled optical depth along the line 
of sight; panel (c) - a cut from panel (b) for 7 = (black solid 
line); the red dashed line corresponds to the dust disc model with 
a simple exponential density distribution law ([T]l; panel (d) - the 
residual map of the observed and the modelled optical depth. 



the data points to vary randomly within their individual uncer- 
tainties. For each model pixel, 100 different sets of data points 
were created and approximated with the two-component model. 
Examples of the black body fits to the data points for three ran- 
domly selected lines of sight (i.e. for three model pixels) are 
given in Fig. |2] The final temperature maps and lOO/^m inten- 
sity maps of the two dust components are presented in Figs. [3] 
and |4] together with the corresponding Icr-error maps. 

Next, the total optical depth T^ax (applicable to light sources 
behind the dust disc) was calculated for each model pixel accord- 
ing to Eqs.|9]and[T0l Only the cold dust component was used for 
calculating the optical depths, thus the calculated opacities refer 
to the diffuse dust component. Since the mass of the warm dust 
component is very low (see Sect. |4|, its contribution to the opti- 
cal depth would be negligible. The derived Tniax,v(^, Y) map was 
four-folded (see Fig.|5^) and approximated with an axially sym- 
metric distribution ^ as described in Sect. |2.2| (see Fig.jSj^)- The 
residual map, indicating deviations between the axially symmet- 
ric model and the four-folded map is shown in Fig.|5}i. 

Approximating the far-IR maps with the axially symmetric 
distribution ^ does no provide an independent estimate for the 
vertical thickness of the dust disc component. We have used a 
rather conventional assumption that the mean ratio of the thick- 
ness of the stellar disc to the thickness of the dust disc is 1 .8 ±0.6, 
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derived by modellin g a sample of nearby edge-on galaxies by 
IXilouris et af] (Il999h. However, in some galaxies the dust disc 
appears thinner (l Bianchill2007h . thus we have made comparative 
calculations also for the cases of a thinner and a thicker dust 
disc (see Sect.|5]for models with different dust disc thicknesses). 
The scalelength of the dust disc was determined during the ap- 
proximation process of the far-IR emission maps. The structural 
parameters for the exponential approximation of the dust disc 
density distribution n{a) are given in Table |2] The dust disc den- 
sity profile (|7]i in the r-direction is given in Fig.|5];, also showing 
that without a correction function /du,st('"), the simple exponential 
function n{a) provides a very poor approximation. 

3.4. Model fitting 

Prior to starting the actual fitting process, a reasonable number 
of different stellar components to be incorporated in the model 
was specified. 

The two most dominant components of M 3 1 are its bulge 
and disc, detectable even with an amateur telescope. In addi- 
tion, recent studies suggest the presence of a dynamically warm 
but also slowly rotating stellar population of stars with mod- 
erate metallicity beginning at the outskirts of the bulge region. 
Here we designate this population as the halo; note that the 
same region has been re ferred to as the bulge dKalirai et al.l 
2006), the extended disc (Ibata e t al.ir2005 l '20qV). the spheroid 
(Brown etal. 2007), and the outer halo (Durrell et al. 200 ft), 
depending on the viewpoint and the focus of each study. 
Furthermore, recent observational studies of distant fields in the 
direction of the minor axis have provided evidence for the ex- 
istence of stellar populations far beyond the conventional bulge 
and halo regions (Irwin et al. 2005; Ibata et al. 2007). The outer- 
most regions seem to be dominated by a faint, diffuse, very ex- 
tended metal poor halo, possibly embedding much of the whole 
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Table 2. Calculated parameters of the photometrical model. 



7 



Population 


ao 




N 


K 


k 


h 




(kpc) 












Bulge 


0.50±0.01 


1.06+0.03 


3.46+0.11 




7.2086 • 10-* 


896.897 


Disc 


6.7+0.2 


0.08+0.01 


1.90+0.13 




6.4723 • 10-2 


26.7655 


Young disc'"' 


9.5+0.2 


0.02+0.003 


0.35+0.02 


0.67+0.02 


1.3334 


1.23802 


Halo 


6.5+0.4 


0.50+0.02 


5.0+0.2 




4.1625 • 10 


31809.6 


Dust disc 


17.0 


0.014 


1.0 




0.5000 


4.00000 


Population 


Lu 


Lb 


Lv 


Lr 


L, 


Ll 




(lO'OLo) 


(lO'OLo) 


(lO'OLo) 


(lO'OLo) 


(lO'OLo) 


(IO'OLq) 


Bulge 


0.63+0.07 


0.82+0.04 


(0.96+0.05) 


(0.93+0.04) 


0.83+0.09 


4.08+0.22 


Disc 


0.61+0.07 


0.90+0.06 


1.36+0.07 


1.47+0.13 


(1.88+0.27) 


5.40+0.52 


Young disc**' 


0.98+0.11 


1.10+0.08 


0.85+0.08 


0.83+0.13 


(1.19+0.21) 


1.09+0.58 


Halo 


(0.79+0.19) 


(0.82+0.11) 


0.23+0.02 


(1.70+0.28) 


0.72+0.06 


(2.32+0.95) 



Notes: The luminosities are corrected for the intrinsic absorption and the absorption in the Milky Way. The luminosities supported by insuf- 
ficient data are given in brackets. Structural parameters for the young disc are only given for the positive component; see Sect. 12. 1 1 for 
determining the parameters of the negative component. The luminosities of the young disc are its true luminosities (L = L+ -(- L_). 



Local Group. Since it has almost negligible surface brightness 
and total luminosity, we do not model this population as a sepa- 
rate component. 

Also, we made no attempt to include the Giant Stream or any 
other of the recently detected stellar streams in our model sep- 
arately. Their contribution to the surface brightness, either mea- 
sured along the axes of the galaxy or averaged along elliptical 
contours is insignificant. 

Consequently, the fitting process was started with three stel- 
lar components: a bulge, a disc, and a halo. However, we soon 
realised that an additional component is necessary for a better 
description of the blue star-forming ring at about 10 kpc from 
the centre, also coinciding with a peak in the distribution of gas 
clouds and other ob jects typical of reg ions of star formation 
jEinasto et al.lll980l; " lTenies & Haud|[l99ll: Itenies et all 11994 . 
The importance of considering clumpy young populations in 
models of the intrinsic dust extin ction of galaxi es is also stressed 
in radiative transfer models (e.g. :Bianchil l2008'). We inserted the 
young stellar population of M 3 1 in our model as a disc with a 
central hole described by Eq.|2j we refer to this population as the 
'young disc' in the following. 

The model parameters q, ao, L, N, and k for each compo- 
nent were determined by a subsequent least-squares approxi- 
mation process. First, we made a crude estimation of the pop- 
ulation parameters. The purpose of this step is to avoid obvi- 
ously non-physical parameters - relation (|4]i is non-linear and 
fitting a model to the observations is not a straightforward 
procedure. Next, a mathematically correct solution was found 
for each component. Details of the least-squares approximation 
and the general model li ng procedure have been described by 
lEinasto & Haul (Il989h . iTenies et al.1 (11994'). and i Tenies et all 
(1998). Uncertainties of the model due to observational uncer- 
tainties and the coupling of parameters were estimated using 
the partia l second derivatives of the sum of the least-square dif- 
ferences dBevington & Robinsonl20d3h . For each parameter, the 
calculated errors only reflect the goodness of the given fit; they 
do not represent the uniqueness of the model. Models with con- 
siderably diff'erent parameters can be fitted to the observational 
data, yielding somewhat larger deviations. 

4. Results 

In the presented model, the dust disc of M 3 1 has been assumed 
to consist of a cold component and a warm component. The tem- 



perature and intensity distribution maps together with their Icr- 
errors for the warm component are given in Fig. [3] It is seen 
that the temperature of the warm dust is 56-60 K in most cases, 
with a typical uncertainty less than 1 K. The highest tempera- 
tures and intensities can be found in the spiral arms, revealing 
a relation between the warm dust and the star-forming regions. 
Figure |4] shows the temperature and intensity distribution maps 
of the cold dust. Temperature of the cold component is usually 
15-19 K with an uncertainty 0.2-1.5 K. In contrast to the warm 
dust, the temperature of the cold dust rises towards the centre 
of the galaxy, up to 25 K. This can be explained by a stronger 
general radiation field resulting from the higher concentration of 
stars. The intensity distribution of the cold component is much 
smoother than that of the warm component, matching well with 
its more diffiise nature. The integral dust temperatures estimated 
bv lGordon et all (|2006|) . 17 K for the cold component and 59 K 
for the warm component, agree well with the average values de- 
rived in our study. 

We have used two-dimensional far-IR maps of M 3 1 for our 
extinction analysis. For illustration. Fig. |6^ shows the Spitzer 
far-IR measurements, averaged along elliptical isoluminosity 
contours. In Fig.|6]3, the calculated 100 /im emission distribution 
is compared to the corresponding observations by IRAS; the low 
resolution of IRAS prevents from seeing the wavy pattern of the 
spiral arms. Figure|6j3 also shows the contributions to the 100//m 
flux by the warm and the cold dust. Except for the very centre of 
the galaxy, contribution by the warm dust is negligible. 

Using Eq. |9] we estimated the ratio of the integral cold dust 
column density to the warm dust column density to be (7 ± 1 ) ■ 1 0-'. 
Since the space density distribution of dust is proportional to the 
column density, the given value is a direct estimate of the ra- 
tio of the col d dust mass to the warm dust mass. According to 
iGordon et al.l (l2006h . the mass of the cold dust component ex- 
ceeds the mass of the warm component by about 10"^ times, close 
to our estimate. The warm dust gives the highest contribution in 
the inner bulge region and in the spiral arms, where up to 1 % 
and up to 0.1% of the total dust mass can be ascribed to the 
warm dust, respectively. This demonstrates the insignificance of 
the warm dust as an absorber of starlight. On the other hand, its 
contribution has to be estimated and subtracted carefully from 
far-IR maps to avoid an overestimation of the extinction caused 
by the cold dust. 
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Fig. 7. The observed (data points) and the modelled, dust-extincted (thick red solid lines) surface brightness distributions of M 3 1 
together with the contributions by the individual stellar components (dashed lines). Thin solid lines represent the intrinsic, dust-free 
surface brightness distributions. Left panels: UBVRL (from top to bottom) surface brightness along the major axis; right panels: 
UBVRI (from top to bottom) surface brightness along the minor axis. 



The derived maps of the temperature distribution and the 
100 yum flux of the cold component allowed us to calculate the 
total optical depth along each line of sight (Fig.|5j)). Within the 
main body of M 31 the total optical depth is Tmax = 0.8-1.2 and 
drops to about 0.3 at the physical limit of the currently available 
Spitzer far-IR data at about 17 kpc. Beyond that, we have as- 
sumed its smooth decrease to zero in our model. The map of the 
total optical depth was used for estimating the intrinsic extinc- 



tion of M 3 1, but it can also be used for estimating the attenuation 
of background objects. 

According to the model presented in Sect. 12.21 the number 
density distribution of the dust disc is represented by an axi- 
symmetric distribution (|7]l. In the r-direction, the number den- 
sity distribution is determined on the basis of the far-IR intensity 
distribution. The derived number surface density distribution of 
dust (which is proportional to r^ax) for M 31 is given in Fig.|5]:. 
According to the present model, an approximation of the dust 
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Fig. 8. Distribution of the observed (data points) and modelled 
(solid lines for the apparent and dashed lines for the dust- 
corrected) colour indices along the major semi-axis. 
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(solid lines for observed and dashed lines for dust-corrected) V- 
/ colour index along the minor semi-axis. 



disc with a simple exponential law using = 1 in Eq. [T] gives 
(flo)dust = 17 kpc as its effective scalelength. The scaleheight of 
the dust disc is set to be 1.8 times the scaleheight of the stellar 
disc and is thus 0.24 kpc. 

The parameters of each stellar component in the best-fit 
model: the harmonic mean radius oq, the axial ratio q, the struc- 
tural parameters A^, the dimensionless normalising constants h 
and k, the relative size k of the hole in the young disc component. 
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Fig. 10. Dust extinction in M31 for U, B, V, R, I, and L filters, 
calculated from the model. Upper panel: extinction along the 
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and the absorption-corrected t/BVTJ/L-luminosities are given in 
Table |2] together with the uncertainty estimates. The calculated 
surface brightness distributions are in good agreement with the 
observed surface brightness distributions. 

Deviations from the observations can be further lowered by 
increasing the number of components. On the other hand, due to 
the uncertainties of the surface brightness data and the moder- 
ate number of available colours, the component parameters are 
degenerate and the best-fitted solution is not unique: a change 
of some parameters of one component can often be compen- 
sated with a change of certain parameters of another compo- 
nent. Without including additional observational data, a higher 
number of components would decrease the uniqueness of the 
model, probably also decreasing its correspondence to the real 
galaxy. The aim of the present paper was to map the effects 
of dust attenuation, thus metallicity determinations in M 3 1 and 
chemical evolution model calculations have not been used for 
laying further constraints on the component parameters. More 
deta ils on the coup ling of the parameters of M 3 1 can be found 
in T enies et alj(ll99 4). At present, the lack of data leaves the halo 
to be the most undetermined component, especially in the U, B, 
R, and L colours. Similarly, because of a slight inconsistency of 
the observational data close to the centre of the galaxy, the V 
and R colours and the corresponding colour index of the bulge 
cannot be determined well. The / colour surface brightness dis- 
tribution is available only along the minor axis. For this reason 
we cannot clearly distinguish the contributions of the disc and 
the young disc to the / colour and the corresponding luminosi- 
ties are rather uncertain. 

The axial ratio of the bulge in our best-fit model is ^ = 
1.06 + 0.03, referring to a slightly prolate shape of this com- 
ponent. This r esult seems to be cont r adicting ear l ier estimates 
q ^ 0.6-0.8 dde VaucouleursI 119581: iKentI 119891: iTenies et all 
I l994|) . However, in those studies the intrinsic absorption has 
not been considered. Extinction is stronger for a given distance 
along the minor axis than for the same distance along the ma- 
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jor axis in the bulge region, causing a decrease of the apparent 
axial ratio. Without the intrinsic absorption, our model would 
yield q - 0.8 as the axial ratio of the bulge, in agreement with 
previous studies. 

To compare the scalelength of the dust disc with that of the 
stellar disc we fitted a single exponential stellar disc to M31; 
the model parameters for the bulge and halo components were 
kept fixes according to the four-component approximation. The 
single disc has a scalelength (ao)stars - 9.2 kpc, thus (flo)dust - 
1 .8 (a o)stars, which agrees w ith the seven late-type galaxies stud- 
ied bv lXilouris et"ail(ll999l) . 

Figure|2]shows the modelled surface brightness distributions 
along both the major and the minor axes, together with the corre- 
sponding observational data. The contribution of the individual 
components and the intrinsic, dust-corrected surface brightness 
profiles are also shown. Note that while moving from U to L, 
the young disc component becomes less dominant, because the 
luminosity of the star-forming ring is dominated by young blue 
stars. In Figs. [8]and|9l the distributions of the modelled and the 
observed colour indices are shown. 

Figure [To] presents the calculated dust extinction as a func- 
tion of radius for different filters. To some extent, these distri- 
butions are model-dependent, but should be applicable with rea- 
sonable precision for correcting any previous or coming obser- 
vations of the surface brightness distribution of M 3 1 . The largest 
uncertainties are related to the derivation of the optical depth of 
the dust disc from the infrared observations. The coiTesponding 
errors in Fig. [10] would be of the order of 20 per cent (for the 
uncertainties of Av see Fig.fTTI). Dust extinction along the minor 
axis is higher than along the major axis due to projection effects. 
The minor axis data represent the closer side of the galaxy, where 
most of the stars lay behind the dust disc, while along the major 
axis, the number of stars lying in front and behind the dust disc 
is equal. This effect is also illustrated by Fig. [7] 

The total luminosity of M31, corrected for the intrinsic ab- 
sorption is Lb - (3.64 ± 0.15) ■ 1O'"L0, coiTesponding to an ab- 
solute luminosity of Mg - -20.89 + 0.04 mag. The dust disc ab- 
sorbs as much as 20 % (or 0.24 mag) of the total B-flux. Roughly 
half of the B-flux is radiated by the discs (25 % by the old disc 
and 30 % by the young disc), about one quarter (22.5 %) by the 
bulge, and one quarter (22.5 %) by the halo. The integrated dust- 
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corrected intrinsic and visible colour indices are given in Table[3] 

The ratio of the luminosity of the spheroidal structures (the 
bulge and the halo) to the total luminosities is 0.45 in B and 0.35 
in V. 



5. Discussion 

In Fig. [TT] we display a comparison of the derived dust extinc- 
tion profi le to some earher resu lts. Assuming a constant dust to 
gas ratio. IBattaner et al] (Il986h have used the H i distribution to 
estimate visual extinction Av inside M31. Considering the dif- 
ferences in methodology and in the availability of data, their es- 
timate matches well with our extinction profile. iHodge & Led 
(11988 ') have used theU-B versus B-V colour diagrams to de- 
termine reddening inside the disc of the galaxy. If converted into 
extinction, these values remain somewhat higher than our esti- 
mates, especially the 10 kpc measurement, which is off by 0.7 
magnitudes; however, also the errorbars are large. IXu& Heloul 
(1996) have estimated the visual optical depth using the IRAS 
satellite 60 and 100 //m maps, an ultra-violet map. Hi map, and 
a model of heating and cooling of dust grains. Converting their 
optical depth estimates to extinction yields a slightl y higher level 
of ext inction than calculated from our model, but IXu & Heloul 
d 19961) have also shown that the calculated optical depth has a 
large scatter across the ga laxy. The most recent analysis of dust 
inside M31 was done bv iMontalto et al.l (12009!) . In their study, 
the Spitzer data have been processed in a way quite similar to 
ours. In both studies the main contributor to the far-IR emission 
is the diffuse dust. Also, the location of the peak attenuation is 
the same, 10-12 kpc from the centre. According to our calcula- 
tions, the attenuation decreases towards the outer regions faster 
tha n towards the inner re gions. 

iBarmbv et all ( |2000|) and iFan et all ( |2008|) have measured 
the mean reddening of globular clusters of M31 caused by the 
dust of the galaxy to be EiB - V) = 0.16 and E{B - V) = 0.28, 
respectively. The former value (corr espondiiig to A y - 0.5) 
agrees with our estimates, whereas Fan et al J ( 12008!) have de- 
rived a somewhat higher extinction (Ay - 0.87) than suggested 
by our model, but once again, the scatter of these reddening es- 
timates is large. 

In Fig. [12] we compare our results for the dust tempera- 
ture and intensity distributions with distributions of molecular 
gas and young stars. The molecular gas distributi on has been 
measured as the intensity of the CO emission by Nieten et al.l 
(i2006) : the surface density distribution of a sample of stars with 
U — V < -0.9 from lBerkhuiisen & Humphreys! ( 11989!) is used as 
the distribution of young stars. To enable a straightforward com- 
parison, the dust parameters have been measured within a nar- 
row stripe along the major axis, similarly to the measurements 
of the distributions of gas and stars. No significant coiTelation is 
seen between the cold dust temperature and the molecular gas 
or the young stars (Fig. [T2h). The temperature of the cold dust 
seems to be determined by the general radiation field of the old 
stars. However, a substantial amount of the emission by the cold 
diffuse dust can also be powered by non-ionising ultra violet 
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Fig. 12. Distribution of the dust properties derived from present 
model (black solid lines) toget her with the distr ibutions of 
CO emission (green do tted line, Nieten et al. 2006 ) and young 
O stars (red histogram, iBerkhuijsen & Humphreys! 19891) along 
the major axis, (a) - cold dust temperature; (b) - warm dust tem- 
perature; (c) - cold dust emission intensity at 100 yum; (d) - warm 
dust emission intensity at lOOyum. 



photons coming from young stars, as shown by 'Popes cu et alJ 
(2000) for the case of NGC 891. On the other hand, a correla- 
tion exists between the amount (intensity) of the cold dust and 
the distribution of the molecular gas (Fig. fT2b). A difference is 
only seen at the centre of M 3 1 where the intensity of the cold 
dust emission is high while almost no gas is detected. A peak at 
the centre is also seen in the intensity distribution of the warm 
dust (Fig.fTSb). Correlation between the intensity and tempera- 
ture distributions of warm dust and young O stars is moderate 
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Fig. 13. Temperature distribution of M31, derived by fitting 
ellipses to the temperature maps. Dark grey area - effective 
temperature; blue and red errorbars - warm and cold temper- 
at ure componen t . The black solid line is temperature derived 
bv lBianchi et alj (|2000|) for NGC 6946. Temperature distribution 
of NGC 6946 is scaled to M 3 1 by taking different diameters of 
these galaxies into account. 



only (Figs. [T2b and[T2h). The offset can be caused by the ten- 
dency to see CO emission from the trailing edge of the spiral 
arm, while most of the detectable 24 fj.m emission (radiated by 
the warm dust) is coming from the leading edge. 

We have modelled the far-IR spectral energy distribution as a 
superposition of two modified Planck functions along each line 
of sight. However, we do not actually know whether the two dust 
components are physically distinctive or perhaps simply a math- 
ematical description of a single dust disc with higher tempera- 
ture variations. To provide a description for the case of a single 
dust component, we have calculated the distribution of the mean 
intensity-weighted temperature of the dust. In Fig.[T3]the mean 
temperature is given as a function of radius. It is seen that the 
temperature decreases rather smoothly from T ~ 32K at the 
centre to T ~ 20K at 7? ~ 7 kpc and outwards. For illustration. 
Fig. [13] also shows the dust temperature distribution for the Scd 
gala xy NGC 6946, dete rmined using a radiative transfer model 
by Bianchi et all (|2000|) ; for this plot the radii of NGC 6946 have 
been multiplied by a factor of 5.9 according to the ratio of the 
half-light radii of these galaxies. It is seen that although the dust 
temperature in these galaxies varies within the same range, it de- 
creases more steeply in M 3 1 indicating a more vivid bulge-disc 
distinction. 

With the help of the presented model, we can also study more 
general dependencies of dust extinction on the inclination angle, 
the total optical depth and the dust disc thickness of a galaxy. 
In Fig. [14] extinction calculations are presented for three dust 
disc relative thicknesses (the axial ratios q): 0.005, 0.014, and 
0.042, designated as 'thin', 'medium', and 'thick', respectively; 
and for three values of the optical depth t: 0.5Tmax, 1 .Or^ax, and 
2.0T,nax, where T„iax is the optical depth as derived for M 3 1 . The 
other parameters have been set according to the model derived 
for M31. The vertical line in Fig. [T4l refers to the inclination 
of M 3 1 . In the upper panel, the extinction is calculated for the 
whole galaxy, in the middle panel for a pure bulge component, 
and in the lower panel for a pure disc component. 
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Fig. 14. Integrated dust extinction inside M 3 1 in V-flux as a 
function of the inclination of the galaxy for a variety of optical 
depths and dust disc thicknesses, r^ax is the optical depth used 
in our model. Upper panel shows the extinction for the whole 
galaxy, middle panel for a bulge component, and lower panel for 
a disc-like component. The vertical line represents the inclina- 
tion angle of M 3 1 . 



In the case of the above-described model of M31, the to- 
tal extinction would be maximal if the inclination angle of the 
galaxy were approximately 88°. At lower inclination angles (a 
more face-on orientation) the total extinction decreases as a re- 
sult of the decreasing line-of-sight optical depth. For an edge-on 
galaxy (1 - cos(0 = 1), the visible area of the dust disc becomes 
negligible, and the total extinction is low. It is seen from Fig. [14] 
that extinction is lower for the disc component than for the bulge 
and increases more rapidly while moving to higher inclination 
values. The extinction maximum occurs at higher inclinations 
for a disc than for a bulge. Similar dependence o f dust attenua- 
tion on inclination was predicted by Tuffs et a ll (|2004|) and are 
now confirmed observationally bv lDriver et alj ( 20071) using dif- 
ferent models and data. 

Figure [14] also shows that the extinction maximum shifts 
slightly towards the lower inclination angles as the optical depth 
increases; of course, the total extinction increases as well. On the 
other hand, the total extinction is almost insensitive to the dust 
disc thickness at low and intermediate inclination angles. This 
is expected, since the extinction of light emitted outside the dust 
disc only depends on the optical depth of the dust disc and the 
dust disc geometry only affects the extinction of light emitted 
inside the dust disc. In most cases, the line-of-sight thickness of 
the dust disc is smaller than that of the stellar components and 



therefore the total extinction does not change when changing the 
dust disc geometry. 

In cases of very high inclination angles, the dust disc thick- 
ness becomes much more important, especially if the optical 
depth is higher as well. However, akeady for M31 and galax- 
ies at lower inclination angles, the line-of-sight optical depth is 
the dominant factor 

The dust model used in the present paper involves certain 
simplifications. Firstly, dust temperatures have not been deter- 
mined on the basis of a complete radiative transfer model. We 
have assumed that the distribution of the dust temperature along 
each line of sight consist of a cold component and of a warm 
component. The temperature maps have been derived by deter- 
mining the corresponding four parameters (the temperatures and 
the intensities of the two components) according to four observa- 
tional data points for each line of sight. Although mathematically 
correct, the derived temperatures may be uncertain since the 
chi-squared approximation is not overdetermined. Secondly, be- 
cause we do not use a radiative transfer model, the radiation scat- 
tering is also ignored. Thirdly, our model does not consider pos- 
sible small-scale dumpiness of the dust. Although the density 
behaviour of the dust disc does not follow a simple smooth ex- 
ponential decrease, having a more complicated structure thanks 
to the empirically determined correction function /dustC?") (see 
Eq.|2]l, possible clumps of dust smaller than the resolution of the 
Spitzer far-IR imaging may be present. The presence of small 
dust clumps would lead to a lower level of extinction than es- 
timated by our model, since the number of unaffected stars be- 
comes larger 

The vertical thickness of a realistic disc with a spiral pattern 
should be varying, but the thickness of the old and the young 
stellar disc is assumed to decrease monotonically in the present 
model. However, the density distribution of the young stellar 
disc has been modelled as a ring-like structure with a central 
density minimum. Since the young disc is thinner than the old 
disc, their interplay gives a slightly varying total effective thick- 
ness for the composite stellar disc. As a result, the ratio of the 
stellar disc thickness to the dust disc thickness is also variable, 
being larger in the regions dominated by the old disc and smaller 
in the other regions. A more detailed study of the structure of the 
spiral arms remains beyond the scope of the present paper. 

The intrinsic luminosity and colour distribution can also 
be used to derive the visible mass distribut ion of galaxies. 
Supplemented with a dynamical model (e.g. Tempel & TeniesI 
2006) and kinematical observations, the visible mass distribu- 
tion can be used for deriving the distribution of dark matter in 
galaxies. We will take a closer look at the distribution of visible 
and dark matter in M 3 1 in a forthcoming paper. 
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